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Fine asymptotic behavior in eigenvalues 
of random normal matrices: Ellipse Case 


Seung-Yeop Lee* and Roman Riser'^ 


Abstract 

We consider the random normal matrices with quadratic external potentials where the associated 
orthogonal polynomials are Hermite polynomials and the limiting support (called droplet) of the 
eigenvalues is an ellipse. We calculate the density of the eigenvalues near the boundary of the droplet 
up to the second subleading corrections and express the subleading corrections in terms of the curvature 
of the droplet boundary. From this result we additionally get the expected number of eigenvalues 
outside the droplet. We also obtain the asymptotics of the kernel and found that, in the bulk, the 
correction term is exponentially small. This leads to the vanishing of certain Cauchy transform of the 
orthogonal polynomial in the bulk of the droplet up to an exponentially small error. 
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1 Introduction and results 


Consider the set of n point particles, {zi,--- ,Zn}, in the complex plane C that interact via two- 
dimensional (i.e. logarithmic) Coulomb repulsion and are subject to the quadratic conhning poten¬ 
tial 

V{z) = \z\^ — tB^e{z^)^ 0 < t < 1. (1) 

The total (electrostatic) energy of the system is given by 


2 1 

= E ^ E ‘“S ^73^. 

id.i i<k ' 1 ^*1 

We consider the canonical ensemble of the particle system, i.e., we assign the Gibbs measure, 

n - n 

Probn({^di=i) n ^ exp ( - iVL;„({zd”=i)) H 

i=i r=i 

= T ( n 1^, - n.,) 

^ j=l 


( 2 ) 


on the configuration space of the n particles. Above, A is a positive constant, dA(z) is the two-dimensional 
area measure, and Zn is the normalization constant. 

The corresponding random process is a special case of Coulomb gas or (3-ensemble (see, for instance, [16]). 
It also describes the eigenvalues of a random normal matrix ensemble [14, 19, 28]. The case of t = 0 is 
called Ginibre ensemble and has been studied in [18, 17, 26]. The case of non-zero t has been studied, 
for instance in [5]. 

The (averaged) density, pn, of the particles, is given by 


Pn{z) :=e(^'^6{z - Zj)\ (3) 

^ i=i ^ 

where the delta function is with respect to the area measure, and the average is taken with respect to 
the Gibbs measure (2). Note that we normalize with 1/A such that the total mass of pn is n/A instead 
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Figure 1: pn for t = 0.25 and n = N = 50 


of one. In fact, throughout this paper, we consider the scaling limit where the total mass of pn is fixed, 
i.e., for a fixed T > 0, N grows with n such that 


T = 


n 

N' 


In Figure 1, one observes that the density is supported approximately on an elliptical shape. In fact, as 
n grows to infinity, the density pn converges to (1/vr times) the characteristic function over 


K =lz 


(Re z)'^ + \ (Ini z)^ <T \ CL 


1 +1 


1-t 


( 4 ) 


Remark. The weak convergence of the density is known [2, 14] for more general class of potentials where 
the support, K, of the limiting density can be an arbitrary compact set in C depending on the potential. 
In [2], K is called a droplet] we will use this terminology throughout the paper. 

Figure 1 shows that the density changes drastically near the boundary of K. Zooming at the boundary 
with a scaled real coordinate ^ normal to the boundary, certain universal shape arises in the density 
profile, see Figure 2. The dashed line in the figure is given by [17] 

where 

2 r°° 2 

erfc(z) = —j= / e”"' dtc. 

Jz 

In [3] the above limit is proven for other classes of potentials. In the next theorem, we show the more 
detailed behavior of the density for the potential given by (1), including the subleading corrections in 
1/A^-expansions. 
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Figure 2: 



Theorem 1.1. Given T > 0, let the compact set K G C be given by (4). Let zq be a point on the 
boundary of K. Let k = k{zo) be the (positive) curvature of dK at the point zq, and dgK (and dgn) be 
the derivative (and second derivative respectively) of the eurvature with respect to the counter-clockwise 
arelength parameter s of the curve dK. Let n = n(zo) be the outward unit normal vector from K at zq 
and a real coordinate. 

Let 0 < n < 1/6. As the positive integer n and the real number N go to infinity while keeping n/N = T 
fixed, we have the asymptotic expansion 


Pn 




= ^ erfc(\/2^) + 
ZTT 

1 f 2 


_^(£2 

^/N 3V27r^ 

- 8£3 + 3£ 
18 


-l)e-2«^ 

, (jdsKf 

V 9k2 



+ O (^iV-3/2+91/^ ^ 


(5) 


where the error bound is uniform over |£| < iV^ and over zq G dK. 

In Appendix A, we show a few numerical plots regarding the subleading corrections of pn- For explicit 
expressions of k, dgK, and d^K, see (88), (89), and (90). 

Remark. The reader may wonder why we choose to write the expansion (5) in terms of the curvature. 
There are two reasons: “simplicity” and “possible universality”. In fact, we conjecture that the term of 
the order 1/y/N in (5) is universally true for a general droplet with a real analytic boundary. The term 
of order 1/N does not seem to hold in general and, therefore, is chosen purely for the simplicity of the 
expression. 


Theorem 1.2. Let 0 < r < 1/2. Given T > 0 there exists constants Cq > 0 and Nq > 0 such that for 
all N > Nq and n = NT, we have 


Pn{z) 


1 

TT 




\pn{z)\ < 


4 


when z G Int(iF) \ Ut^n-, 
when z G Ext(iF) \ Lfr,N- 














We denote the neighborhood of dK of width 2N by 

Ut,n = {^0 + XTa.{zo) I zq G dK, |X| < = {z \ dist(z, dK) < 


Choosing 0<r<z^<l/6, Theorem 1.1 and Theorem 1.2 together give a uniform estimate of pn over 
the whole complex plane. 

It is also interesting to look at the kernel function defined by (see, for instance, [16]) 

n—1 

Kniz,w) = Y, ( 6 ) 

j=0 

where the polynomial pj of degree j is the normalized orthogonal polynomial satisfying 


[ Pjiz)pkiz)e = 6jk. 

Jc 


(7) 


The following theorem describes the asymptotic behavior of Kn- 

Theorem 1.3. Given T > 0, let the compact set K C C be given by (4). Let Vk be a compact subset of 
Int K and Uk a neighborhood of K. Then there exist e > 0 and 5 > 0 such that, as n = NT grows to 
infinity, the following asymptotic expression holds uniformly over all w and z as specified below. 


1 

~N 


Kn{z,w) 


— exp (—y|z — exp ^zA^Im^ 2 ;r(;——(l + C>(e w ^ Vk, z ^ Bs{w), 

' 0{e-^^), weVK,z^Bsiw), 

.0(e-^^), w^Uk,z€C, 


The notation Bs{w) stands for the disk of radius 6 centered at w. We note that the above holds when z 
and w are exchanged. 


Let us consider more general potential given by 

V{z) = \z\^-Q{z)-W)- ( 8 ) 

where Q : C —)■ C is an analytic function (in a neighborhood of K) such that lim|^|^oo V(z) / log \z\ = + 00 . 
Then the first asymptotic formula in the above theorem can be written more generally as (9) below. For 
a general potential the leading behavior of the kernel is known in [2, 6]. Our new observation is that, 
for Q{z) = tz^l2, the error bound in (9) is exponentially small in N. We conjecture that this is true for 
more general Q. In the next theorem, we connect the behavior of the kernel with the behavior of certain 
Cauchy transform. 

Theorem 1.4. Let V be given by (8) and let Q be analytic in a neighborhood of z € C. Let pn and Kn 
be the corresponding orthonormal polynomial and the kernel as defined in (7) and (6) respectively. In the 
scaling limit where n = NT goes to 00 for a fixed T > 0, let Kn satisfy the asymptotic behavior: 

^Kniz,w) = ie-fl"-H%*^i“("--QW-^)(l + 0(e-^^)), (9) 

N TT ^ ^ 
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uniformly for w £ Bs{z) for some e > 0 and 6 > 0. Then there exists e' > 0 such that 


\Pn{w)\^e 


2^-NV{w) 


Z — W 


-d74(t(;) = 0{e 


-Ne'\ 


(10) 


Corollary 1.5. When Q{z) = tz^I2, the Cauchy transform at (10) vanishes exponentially fast in N 
when z £ Int K. 


Though it follows directly from Theorem 7.7.3 in [2] that the above Cauchy transform vanishes in the 
leading order, it is new (at least to us) that the Cauchy transform vanishes at all orders in the 1/N 
expansion. The Cauchy transform in (10) is related to the Cauchy transform that appears in the Dbar 
problem [21] through the following identity. 


1 

Pn{z) 


\pn{w)\^e 


2 „-NViw) 


Z — W 


-d^(t(;) 


Pn{w)e 
z — w 


d^(rc), 


which can be verified by the following equation that is immediate from the orthogonality of PnS. 


r Pnjz) -Pnjw) 

Jc z:-w 


Pn{w)e-^^^^^dA{w) = 0. 


In the remainder of this paper, we will present the proofs of the four theorems. In Section 2 we list the 
notations and a few useful (known) facts. In Section 3 we give the asymptotic behavior of the orthogonal 
polynomial Pn{z) using the properties of the Hermite polynomials. In Section 4 we prove Theorem 1.3. 
In Section 5 we prove Theorem 1.1 and Theorem 1.2. In Section 6 we prove Theorem 1.4. In Section 7 
we discuss several issues including the expected number of particles outside the ellipse. 


2 Useful definitions and facts 


The foci of the ellipse dK lie at ±To where 


Fn = 


AtT 


1-F' 

For 0 < t < 1 we define the conformal map (D is the unit disk) 

4> : C\ closD -£■ C\K by (f){u) := H—^ . 

One can check that (/>(5D) = dK. 

The inverse function gives the conformal map 

iP :C\K ^C\closB by V’(-z) = - F^/z^y 


( 11 ) 


( 12 ) 
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The branch of the square root is taken such that = 2^/1 zjFQ + 0(1) as z goes to oo. It will be 
useful to denote the logarithmic capacity of K by 


cap(it') = lim ^ ^ = \ - -(13) 

^ 2^00 ^{z) 2\/t V 1 -^ ^ 

The Joukowsky map (j) is in fact a conformal map between the larger domains, i.e. 

</.:C\closSy^(0)^C\[-Fo,To], 

where Br{x) stands for the disk of radius r centered at x G C. Therefore, ijj is analytically continued to 
the conformal map 

V^:C\[-Fo,Fo]^C\clos5^(0). 

Throughout the paper it is convenient to consider the expression 

to be analytic on C \ [—Tq, Fq] such that the value at oo is one. 

The Schwarz function S for an analytic curve T in the complex plane is defined to be the unique analytic 

function in a neighborhood of T which fulfills S{z) = z on T. More details about the Schwarz function 

can be found, for instance, in [8]. 

For the ellipse, dK, the Schwarz function S : C \ [—Fo,Fo] —)■ C is given by 

S’(z) = fz + ^^ ^1 - Y^l-F2 /z2^ . (14) 

This can be derived by the following identity 

S{z) = (j){l/^p{z)), 

which holds for a general simply-connected droplet. It is simple to check that S'(z), as defined by the 
above equation, is analytic and satisfies S{z) = z on z £ dK. 

Remark. For t = 0 we have (j){u) = ^/Tu, ip{z) = zj\fT^ cap(iF) = \fT and S'(z) = Tjz. 

Let pj be the orthonormal polynomial of degree j satisfying 

[ Pj{z)pk{z)e~^^‘'^'^dA{z) = 6jk, lim > 0. 

The following statement is known [29] (see also [27]). 

Theorem 2.1. Let 0 < t < 1. The orthonormal polynomials {pj}j=o,i, 2 ,- are given by 

Pjiz) = 


\ i/2 

j (l-t2)V4 



(15) 
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where Hj is the Hermite polynomial of degree j given by the generating function 


^2xz—z‘^ 


The leading coefficient of pj is given explicitly by 




j=0 


X] — . 

J'- 


7 ,' := lim 


Pj{z) _ 


■j+1/2 


z^oo 


V^- 


(16) 


We note that 7 j and pj (and many other related quantities) also depend on the parameter N. 
Remark. If t = 0 the orthonormal polynomials are simply the monomials where 7 j = . 

Lemma 2.2. Given a fixed r € Z, as n = NT goes to oo for a fixed T, the following expansion holds. 

1 + 6r + 6r^ 


1n+r 


WJLi 


1/4 


o-N^ 


1 - 


V27rV y (cap(ii:))2'-+i V' 24rw 

where cap(K) is defined at (13), and the Robin constant I is defined by 

£:=r(21ogcap(is:)-l). 




(17) 


One can prove this lemma by direct manipulation of (16) using Stirling’s formula for the factorial. The 
leading asymptotic expression in the lemma also appeared in [4] which is believed to be universal. 


Let us define the pre-kernel JCn by 


K,n{z,w) = Kn{z,w)eyrp ( ) exp 


N. 


(^N lTa.(^ — zw + ^z^ + ^ 


X] Pi~ 

j=o ^ ^ 


to i' _9 

— i\ \ zw -z- w 

2 2 


(18) 


Proposition 2.3. Let n G N and w,z G C. Then the pre-kernel K,n fulfills the following Christoffel-Darboux 
type identity. 

^ ^^”(z,iu) = --^L={tpn{z)pn-i{w) - Pn-i{z)pn{w)) eyrp l^z^ - (19) 


N dz 


This identity is known [7] in more general cases in the context of the two matrix model and bi-orthogonal 
polynomials. Below we reproduce the proof tailored to our case (See also [27]). 

Proof. We set p-i{z) = 0. From Theorem 2.1 we have (see, for instance, [27, 13]), for A: G N U {0}, 

zpkiz) = rk+iPk-eiiz) + trkPk-iiz), 

1 ( 20 ) 
j^dzPkiz) = trk+ipk+i{z) + rkPk-i{z), 















where we set and denote Pk{z) = Pk{z) The second relation can be derived from 

p'k{z) = yjNk{l - t‘^)pk-i{z). 

Using the definition of the pre-kernel (18) and the recurrence relations (20), we get 


^Nzw 


N 


/ \ ^Nzw ^ _ 

dz[K:niz,w)j = -J^dz ( J^Pj(z)pj(u;)e 


—Nzw 


n—1 / 

^ f {trj+ipj+i{z) + rjpj_i{z))pj{w) - pj{z){rj+ipj+i{w) + trjpj_i{w)) 
j=0 ^ 

n—1 n n—2 

{trj+ipj+i{z) + rjpj_i{z))pj{w) - 
j=0 j=l j=0 


= trnPn{z)pn-l{w) - T nPn-l{z)pn{w). 

This proves the proposition because 

Using Kn or /C„, the average density pn in (3) is given by (note that it is normalized to T) 

Pniz) = ^Kniz,z) = ^ICniz,z). 


Corollary 2.4. Let n be a positive integer. The density pn satisfies 

dpn{z) 


dz 


= {tpn{z)pn-l{z) - Pn-l{z)pn{z)) 


Vi — 


□ 


( 21 ) 


( 22 ) 


3 Asymptotic Expansion of Orthogonal Polynomials 


The orthonormal polynomials, {pk}, are rescaled Hermite polynomials, {Hk}, where the argument of 
Hk is multiplied by the constant '\/N{l — t‘^)/{2t), see Theorem 2.1. The leading asymptotics for the 
Hermite polynomial has been proven by Plancherel and Rotach in [24]. See also [9, 10, 12, 30]. 

We will need the asymptotics of pn to get the asymptotic behavior of the kernel and the density pn using 
Proposition 2.3 and Corollary 2.4. Especially, we will use the (subleading) behaviors of to get the 
(subleading) behaviors of the kernel and the density. 


The following lemma is shown in Appendix B. An alternative, self-contained method to calculate the 
asymptotics of classical orthogonal polynomials can be found in [23]. 

Lemma 3.1. Let n = NT and r G Z. As N goes to infinity, the following asymptotie expansion holds 
uniformly over any compact subset ofC\ [—Fo,Eb]. 


Pn+r{z) 



1/4 ^_ 

fi{zY y/fifiz) exp 


ng{z) - ^ + j^hr{z) + 0 



(23) 
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where if: is given in (12), i in (17) 


and 


g{z) = log z +log + + 

h (.\- - + _ 1 + 6r + 6r^ _ 5^ _ 

~ T\8{z^ - F^) 24v'l - F^/z^ A8{z^ - F^)^l - F^/zV ' 


(24) 

(25) 


Though the function g (24) has the logarithmic branch cut due to the term log z, the expression exp{ng{z)) 
is analytic in C \ [—To, Fq], The expression of hr (25) is analytic in C \ [—To, To]. 

Remark. If t = 0 we simply have g{z) = logz and hr{z) = 

The leading term of the asymptotic formula (23) for the orthogonal polynomial has been conjectured in 
[1] and appeared in [4] for different K. For general K, g is given by the complex logarithmic potential 
generated by the uniform measure on K. 

We will also need the following, rather crude, asymptotic behavior of 


Proposition 3.2. Let n 
whole complex plane. 


NT grow to oo for a fixed T. 


The following estimate holds uniformly over the 


The same bound holds for pn-i. 


Proof. One can use the uniform asymptotics stated at Theorem 1.34, Theorem 1.49 and Theorem 1.60 in 
[11]. These theorems are written for monic polynomials (denoted by Trfis in the cited paper). Therefore 
we need the behavior of the normalization constant 'yj in (16) (which is different from 'jj in the cited 
paper; our normalization is with respect to the area integral) which is given in Lemma 2.2. 

The error bound comes from the behavior of pn near the focal points, see the equation (1.58) 

in [11], and from the behavior of 7 ^ in Lemma 2.2. The former contributes by and the latter 

contributes by □ 


Lemma 3.3. As N goes to 00 while n = NT, there exists e > 0 such that 

p„( 0 ) = - + O(e-'^^). 

TT 

Proof. At the origin the Hermite polynomials give (see for example [20]). 

, J 0 if m is odd, 

\ (—l)”*/^m!/(^m)! if m is even. 

Using this in (15) and (21), the density at the origin is given by 


1 /l- +2 ^ 2 J 

Pn(0) = ^ E \Pii0t = 5: 

j=0 1=0 


t\^^ {2iy. 


(/!)2- 
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Note that the power series in the right-hand side is a truncated Taylor series from 


E 




{2l)\ 


(/!)2- 


(26) 


Choosing e > 0 such that x = < 1, we have 


^Ne 


Pn{0) 

vr 


TT 


^ ^ne/T ~ 

VT^ 


E 


«=LVJ+i 


(21)1 


il\y 


< 


IT 


E 




XN2i (20! 

2/ Jlif 


< 


Vi^ 




where, in the penultimate inequality, we used the index I in the summation is > n/2 and therefore 

gne/T < g2k/T^ □ 


4 Kernel in the bulk 


We define Q : C 


by 


n{z) = Viz)-2TRegiz) + i. 

In fact, the value of i (17) has been chosen such that 17 = 0 on dK. One can show that, for any zq G dK, 


0(z) = |z|2-|zop-2Re [ S{C)dC, 

J ZQ 


(27) 


by taking the derivative and using the explicit forms of g (24) and S (14). (The independence of 11 on 
zq can be seen from 5(zo) = zq which is obtained by taking the derivative of H.) Since ^(z) has a purely 
imaginary jump on [—T(),Tb], H can be defined everywhere, with no discontinuity on [—To,Tb]. The 
following lemma addresses the global minima of 0 (see [27] for another proof). 

Lemma 4.1. For z € C we have 0(z) > 0 and the equality holds only when z G dK. 

Proof. Let us define 11 : C —^ M by 

/ U 

(f{l/v) d^(7;). 

Since cj) : C\ (closR^(O)) —)> C \ [—To, Fq] is conformal, it is enough to show that 

{u G C \ 5^(0) I ll(n) < 0} = 5D. (28) 

Comparing the functions g (24) and V’ (12) we see that 

t . / 2\/l 


g{Hu)) = 




- log 


'tpiz)Fo 


. z=(l){u) 


t , f2Vi 
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A straightforward calculation gives 


Q{Re^(>) = -R'^ + {R^ - l)tcos(20)) - 2riogi2. 

The partial derivatives in 9 and R give 


(9f2(i?e*®) 

¥e 

5^ 


2T 

—^(i? - R~^ f‘t sm.{2e), 

2T(ir^-i;-3) ^^2^^2_ ^. 


(29) 

(30) 

(31) 


The hrst expression vanishes only when either i? = 1 or sin(20) = 0. When i? = 1, the expression (31) 
vanishes. One can also check that 11 = 0 when i? = 1. When sin(20) = 0, the expression (31) vanishes 
only when R = y/i and R= 1. 


We will show that 11 > 0 for R = y/t. There we find 


r 2r(l — l)(cos0)^ 2T(1 — 1) 

HO/te*® =-^- R - Tlogl >-^- - - Tlogl, 

1+1 1+1 


0 < 1 < 1 . 


The right-hand side goes to zero as 1 goes to 1, and it is strictly positive for 0 < 1 < 1 because its 
derivative is strictly negative, i.e. 


r(l -1)2 
1 ( 1 + 1)2 


< 0 , 


0 < 1 < 1. 


Since 11 goes to +oo as R goes to oo, we proved (28). □ 

Lemma 4.2. Let Vk be a eompact subset oflntK and Uk cl neighborhood of K. Consider the limit where 
N goes to oo while n = NT. Then there exists e > 0 such that the following asymptotic behavior holds 
uniformly over all z as specified below. 


Pn{z) 


- + 0(e-^^), 

TT 


z e Vk, 

z G Uk- 


(32) 

(33) 


Proof. From Corollary 2.4 and Proposition 3.2, there exists a constant C > 0 such that 


dpniC) 

df 


< 


(|Pn(0Pn-l(0| + i|Pn-l(0Pn(0|) 


^g-WRey(0 

y/fvru 




(34) 


Since H is strictly positive over a compact subset of Int K by Lemma 4.1, the right-hand side is uniformly 
bounded by where e > 0 can be chosen to be smaller than the minimal value of H over the 

compact subset. To prove (32) we only need to integrate d^pn from the origin to the point z via straight 
path and use the behavior of Pn{0) in Lemma 3.3. 


Using the same idea, the estimate (34) also holds over the complement of a neighborhood Uk of K. Again, 
an integration starting from a fixed point, say zi, in the complement of Uk gives that 


Pn{z) = c{N) + 0(e ^^) 
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where the error term is uniform over the complement of Uk and c{N) is given by Pn{zi)- We note 
that, even though the complement of Uk is unbounded, the integration converges because grows 
quadratically as |^| —)• oo. Then we must have c{N) = 0(e~^'^) because, otherwise, PndA diverges. 

This proves (33). □ 


Now we prove Theorem 1.3. 

For any 2 ; and w the pre-kernel is given by 

^JCn{z,w) = ^JCn{w,w) + ^ J 

The upper bound of the integrand comes from Proposition 2.3 and Proposition 3.2: 


1 

N 


die. 






< 


Vf 


< c 


Vt 


(t \pn{OPn-l{w)\ + \pn-l{C)Pn{w)\J exp - NRe(^^W - 


,Ar5/6 


[“ NRe(^^w - - Tg{i) - Tg{w) + 

for some C > 0 and for all N sufficiently large. The exponent can be further simplified by 
Re(^w - - Tg{^) - Tg{w) + dj =- w\^ + ^^{0 + 

By Proposition 4.1 we can define e such that 


(35) 


(36) 


(37) 


0 < 2e < min Q,{w). 
wGVk 

When ^ = w £ Vk, the right-hand side of (37) is greater than 2e. Then there exists sufficiently small 6 
such that when |^ — t(;| <5 the expression (37) remains bigger than e for any w £ Vk- This implies that 
the upper bound (36) is given by uniformly over w £ Vk- The first case of Theorem 1.3 follows 

from (35) and Lemma 4.2. 


For the second case oi z ^ Bsiw), we will use 


1 

N 


dlCr. 






which can be noticed from (37). 
straight path from w to z, 

l/C„(z,u;)e-fl^-H^ 


exp (-fie - e G C, (38) 

Using this estimate, we get, by integration of d^JCni^,w) along the 

= ^)Cn{w, u;)e-f 1^-"’!^ + z ^ Bs{w)- 


We can use Lemma 4.2 to evaluate ^)Cn{w,w) for w £ Vk- Going back to the original kernel Kn{z,w) 
by (18) proves the second case of Theorem 1.3. Again, as in the proof of Lemma 4.2, even though 2 ; can 
be in the unbounded C, the error bound above is uniform since the error in (38) has U(^) in the exponent 
which grows quadratically for large |^|. 


The last case when w ^ Uk follows from choosing e such that 

0 < 2e < min Q,{w). 
w^Uk 

The rest of the proof remains the same as in the second case. 
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5 Density near the boundary 


We will prove Theorem 1.1 using the Christoffel-Darboux formula in Proposition 2.3 and the asymptotics 
of the orthonormal polynomials Pn obtained in Section 3. 

For any complex number v = Vx + ivy let us denote 

d^fiz) = ivd + vd)fiz) 


where d = 4- = U 4- - iS\ and d = 4^ = U 4- +i4r.]- 


a 


dz ~ 2 \dx ‘'dy 


a _ 1 a 


a 


az 2 \ ax dy 


Let 7 : M —)• dK be the counter-clockwise parametrization given by 

7(0) = 

The (outward) normal vector, n, and the tangential vector, t, at cj){e^^) are given by 

7'(0) 


t = m = 


iy(0)i' 


We will denote the derivatives in normal direction and tangential direction by and dt respectively. 

We denote the derivative with respeet to the arclength by dg, i.e. dg = Obviously dg = dt- 

Generally, 9^ / 5^ for /c > 1. 

Proposition 5.1. Let z G C \ [—Fo,Fo] o,nd v G C. As n = NT goes to oo, we have the asymptotic 
expansion 


»vP„W=(y)‘''e-™''>Re 




\/r(vt - v) \'il^'{z)\ ( h^{z) + h-i{z) 


\/i -y i){z) 


N 


+ Ar2 


(39) 


The error bound holds uniformly on any compact subset o/C \ [—Fo,Fo]. 
Proof. From the Christoffel-Darboux formula (22) we get 

VTe-NViz) 


dvPn{z) = 2Re(^(vt - \r)pn{z) Pn-l{z)j ■ 


(40) 


From Lemma 3.1 we have 




and 






where the error terms hold uniformly on any compact subset of C \ [—Fo,Fo]. Putting these into (40) 
proves the proposition. □ 
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5.1 Expansion of 12 near the boundary 

Recall that dg is the derivative with respect to the arclength and that n is the (outward) normal vector 
at the boundary dK. The (signed) curvature of the curve dK is given by 

d % 

K = — arg n = —idg log n =-(41) 

ds n 

On the boundary we find the following relation between n, k and S. 

Lemma 5.2. The following identity holds for zq G dK. 

n^Sfzo) = —1, n^5''(zo) = 2 k, ii^S”'{zo) = —— lidgK. (42) 


Analogous relations can be found in [ 8 ]. 

Proof. Let 7 : 0 i—)• </>(e®®) be a parametrization of dK. Then we get the tangential vector t by 


t = 


del{G) iu(j)'{u) 


\y{e)\ |(/>'(u)| 


= #'(^o)|^^, zo = </>(e*®). 


The Schwarz function satisfies S{(j){u)) = (j){l/u) in a neighborhood of dK. Taking the derivative in u 
we get 


S {(p{u))(j) [u) = -— 




and, when zq G dK, we have 


S'izo) = -- 




ilj'{zo){if{zo)Y n(zo)2 
using n = —it. This proves the first equation of (42). 

Taking the (9<j-derivative on both sides of (43) we get 


(43) 


iii{zQ)S"{zo) = 2 


n( 2 ;o)' 


rin(zo)K 


(44) 


where we used dg = ind — ind on the left-hand side, and (41) on the right-hand side. This proves the 
second equation of (42). 

The last equation is obtained by taking ds-derivative of the second equation in the lemma. (When 
repeated, such process can yield the relations with higher order derivatives than are shown in the lemma.) 


3 n^(in)KS''' -|- = 2 dgK. 

Using the second equation in the lemma to replace S", one gets the last equation. 


□ 


Remark. Lemma 5.2 and the following lemma are both true for more general droplets. Note that in the 
proof we never use the fact that K is an ellipse. 
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Lemma 5.3. Let zq G dK. Let X,Y G R be defined by 

Xn.{zo) + iYn{zo) = z - zq- 

For any n such that 0 < u < 1/G, as N grows to infinity the following asymptotic expansions holds 
^-Nniz) ^ ^-2NX^ /1 ^ ‘^{X^ - 3XY^) -^{X^- 6X2y2 + y4) 

Y 3 2 

-{X^ - 3XY^f + _ 




3 ' '9 

uniformly over \z — zq\ < and over zq G dK. 

Proof. Using = n5 + nd and dt = ind — ind and the first equation in (42), one finds that 

dnLl{zo) = dtLl{zo) = dndtil{zo) = d^n{zo) = 0 . 

As an example, let us prove the last equality. 

d^Q{zo) = {ind — i\id)'^Q.{z)\^^^^ = (incf — irid) (^rvifz — S{z) — i\i{z — 


(45) 


(46) 


Z=Z0 


= (n25'(z) + n25'(z) + 2) 


= 0 , 


z=zo 


where, in the last equality, we used the first equation in Lemma 5.2. Similarly, we get 


d'^Vl{zo) = (^-n^S"(z) - vfiS'{z) + 2^ 


= 4, 


Z = ZQ 


(47) 


using the first equation in Lemma 5.2. 

In the derivatives of order three or more, the gaussian term \z\'^ in n(z) does not contribute. So we have 

5vi5v2 • ■■dy;^^LL{z) = -25vi5v 2 • • •^v^.Re (^j S{C)dC^ = -2Re (xiV 2 ■ ■ ■ (z)^ , 

for k > 3. Using the second and the third equations in Lemma 5.2, one can obtain the following identities 
for the third order derivatives: 


(48) 


dln{zo) = -dnd^Q{zo) = -4:K, d^Qizo) = -dldtLl{zo) = 0, 
and for the fourth order derivatives: 

d^n{zo) = d^n{zo) = -dld^Q{zo) = 12k^, 9i^9tO(zo) = -dnd^Ll{zo) = -AdgK. (49) 


Since both the real and the imaginary parts of 12 are real analytic, we have 

^ ^ {Xd^ + Ydt)>^n{zo) ^ 


k=0 


k\ 


(50) 
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uniformly over \z — zo\ < S and over zq G dK, for a sufficiently small 6. Putting (46)(47)(48)(49) into 
the Taylor expansion, we get 

n{z) = 2X^ - —- 3Xy2) + '^(^4 - 6 X 2 y 2 ^ Y^) - - XY^) + 0{\z - zo\^). (51) 

3 2 3 

Multiplying the above by N and exponentiating, we get 

exp ( - Nn{z)) = (1 + N—{X^ - SXT^) + N^ — {X^ - 3XY^f 

y 3 9 

- N—{X^ - 6X2y2 ^ Y^) + N^^{X^Y - XY^) 

2 3 

+ O {N\z - zolt N^\z - ^o|^ N^\z - ), 

uniformly over \z — zo\ < and over zq G dK, for a sufficiently large N because the exponential 

function converges uniformly in a bounded region. The lemma follows because, among the error terms, 
N^\z — zo\^ contributes the most. □ 

5.2 Proof of Theorem 1.2 

Theorem 1.2 that we prove in this section is a stronger version of Lemma 4.2. The theorem gives an 

estimate in a neighborhood of dK that scales arbitrary close to the boundary as N goes to inhnity. 

Proof. From (51) we have 

n{z) = 2X‘^ + 0{\z-zof), (52) 

uniformly over \z — zo\ < 6 and over zq G dK, for a sufficiently small 5 > 0. 

We denote the neighborhood of dK of width 26 by 

Us = {zo + Xn{zo) I zo G dK, |X| < 5}. 

For any z € Us\ Ur,N there exists zq G dK such that z = zo + Xn with < |X| < 6. From (52), 

we have P(z) > on z G 14 \ Ut-,n for any sufficiently large N . This 

leads to 

^-Nn{z) ^ g-7V2- q^z£Us\ Ur,N- (53) 

As in Lemma 4.2 we now use (34) to get a uniform bound for the derivative of the density on Us \ Ut^n- 
We get the density at the point z = zq + Xn ^ Us \ Ur,N with zq G dK and < |X| < <5 by 

integrations from starting points zq — 6n € IntiL or from zq + (5n G ExtiF. When N is large enough, 
from (34) and (53), the contribution from these integrations are bounded by 

\Pn{z) — Pn{zo =t (^n)l < 6 ■ when z £ Us \ V-.Af- 

From Lemma 4.2 we know that at the starting points zq^ 6n £ dUs the density equals l/vr + 0(e~^^) 
on Int(iF) \ Us and 0(e~^'^) on Ext(A') \ Us when e > 0 is choosen small enough. □ 
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5.3 Proof of Theorem 1.1 


The proof of the following three lemmas are in Appendix C. 

Lemma 5.4. Let zq G dK. Let X,Y G'R he defined by 

Xn(zo) + iYn{zo) = z - zq- 
The following expansion holds as z goes to zq 

(l + 0i,o(^o)^ + 0o,i(^o)T++ 
fi{z) 'fi{zo) V 2 2 

+ 0{\z- zop) , 

uniformly over \z — zo\ <6 and over zq G dK, for a sufficiently small 5. We used the definitions: 


0i,o(2;o) = -n, 


^ , 1 dgK , ,,, ,, 

© 0 , 1 ( 20 ) = o-*1^ (^o)l, 

O tv 


© 1 , 1 ( 20 ) = + 2iK\ip'{zo)\, 


6k^ 6k 6k 

^ , N 2, (dgK)'^ 9 d'?K , ,,9 , ,,dsK 

© 0 , 2 ( 20 ) = -- ,^2 “ 


9 

Lemma 5.5. Let zq G dK. Then 


Re(/io(zo) + /i-i( 2 o)) = ^ — 

12 24 K 

Im(/io(2;o) - h-iizojf _ dgK 

6k 


IV’'(2o)| 


Lemma 5.6. Let zq G dK. Then 


VT{nt - n) \i:'{zo)\ 
Vl - t‘^ ^^( 20 ) 

VT{m + ii) \'ip'{zo)\ 


= -1 + 


id. 


3|V;'(2o)|k’ 


K 


vi - f^ ^^(20) iV’'( 2 o)r 

Proposition 5.7. Let zq G dK. Let us define X,Y G'S. by 

Xn(zo) + iYn{zo) = z - zq. 

For any 0 < u < 1/6 the following asymptotic expansions hold. 
dr.Pn{z) = (^-2 -k {INX^ - 2X) - k^ {^N^X^ - InX^ + 2X‘^) 


(54) 


+ 


d^K A{dsKf 


6k 


9k2 


v2 1 / K 


N \ 6 


{dsK)‘ 

9k2 




dtPniz) = (^-2 kY + k^ {4NXY^ - ^NX^Y + 4XY) + - Y^'^ 

- - + O ( Ar-3/2+7i.\ ^ 

A^3|i/'(zo)| V J)^ 


(55) 
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where the former is evaluated for Y = 0. The error bounds are uniform over {z \ \z — zo\ < N 
and over zq £ dK. 


Proof. We use Lemma 5.4 and Lemma 5.6 to get 

VT{m-n)\-ip'{z)\ ( hQ{z) + h-i{z) 


Re 


\/l - t^ i’iz) 


1 + 


N 






9k2 


6k 


(56) 


_ Re(feo(zo) + L_i(zo)) dsKlm{ho{zo) - fe-i(zo)) 
N ^ 3k\'iP^{zo)\N 

We use the above with Lemma 5.3 in equation (39) with v = n to get 




dnPniz) = - ^1 + « (|iVX3 -X)+^ (x2 - + iN^X^) + 


2 

9k^ 6k 


+ 


Re(ho(zo) + h_i(zo)) d^K Im(h_i(zo) - ho(zo)) (p/jY 9 ^- 3 / 2 A 
N 3\'tlj'{zo)\K N J 


One can use Lemma 5.5 to finally get (54). 

Similar steps using Lemma 5.4, Lemma 5.6 and Lemma 5.3 in equation (39) with v = t give 

dtpniz) = (^-2 kY + {ANXY^ - ^NX^Y + AXY) + dsK - Y^^ 


2 k : lm{h-i{zo)-hojzo)) / 7 ^_ 3 / 2 ^ 
\f^'{zo)\ N ^ ^ ^ 


(58) 


Again using Lemma 5.5 we get (55). 


Let us remark about the different powers in the error terms in the two cases. For the latter case, dtp, we 
need one less order in the expansion of because the leading term, —2kY = is already 

of less order than the leading term of dnPn- C 


Lemma 5.8. Let zq £ OK. We have, for an arbitrary e > 0, 


1 


Pn{zo) =-- 


2 vr 3\/2tt^N 

uniformly over zq £ dK. 

Proof. Let 6 < v < 1/6. Choose r such that 6 <t < v and let 

z.^ = Zq- 
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for sufficiently large N such that zi G Int K. 


fO 

Pn{zo) - Pnizi) = / dnPnizo + Xn) dX. 

J-N-1/2+T 


Proposition 5.7 gives a uniform estimate of dnPn on the domain of integration. We can replace dnPn by 
(54). Then we only need to perform the integrals of the type 


x4-2^^"dx 


/_Ar-i/2+T 


1 (-ly T. (j+A 
2(2X)0+P/2 \ 2 J 




j G NU {0}. 


The error term is obtained by substituting X = '/NX in the integral and using I’Hopital’s rule, 

XJe-2X^(iX (_iy^Arr(i+i)-ie-2Af"" 

T J—CXD T J—OO T V ^ ^ 

— ITTI - = — nn - = ITTl - 

jv^oc Arfc-(i+i)/2e-2iV2- iv^oo - 4rX2^+fc-i)e-2^"" 

= lim 

4 N^oa 

which converges to a non-zero constant if k = r(j — 1). 

The straightforward calculation of the integration of (54) gives 

Note that the term of order 1/N vanishes. We will discuss a deeper reason for this in Section 7.1. 

The theorem is proved because, according to Theorem 1.2, pnizi) is I/tt up to an error given by 
(9(jV®/®e“'^ ^). The uniformity of the bound over zq £ dK follows from the uniformity of Proposi¬ 
tion 5.7. □ 


We will prove the following slightly stronger statement than Theorem 1.1. 
Theorem 5.9. Let zo G dK and 

z - Zo = n ^ , ^,? 7 GM. 

For any 0 < v < 1/6 the asymptotic expansion 


Pn 



/ + ir])n\ 

Vn ) 


erfc(\/2^) 

2^r 



^ (e-i 

VN\ 3 

\ 







holds uniformly over {z | |z — zo| < X^ ^/2| g dK, where we define 




2 2/ - 8/ -h 3C + 36Cr?2 - + 18f,/ 
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+ 


{dsKfi 


9^2 


d^K\ dsK/ 
12k P 3 V 


2 S 

77 — 77 — 77 


20 






















Proof. The proof is by direct calculation where we integrate dnPn from zq to zq + n^/VfV and use 
Lemma 5.8 for the value of pn{zo). More explicitly we integrate d^Pn using the expression in (54), over 
the integration variable X from 0 to f/y/~N, i.e. 


Pn 




du.Pn{zo + Xn) dX. 


Once we obtain the value of pn at zo+f, n/ \/N, we integrate dtPn from ZQ+f n/\/iV to the final destination 
zo + fn/^/N + rj t/y/N. Again, this is done by integrating the expression in (55) over the integration 
variable Y from 0 to i.e. 


( %/JV ,/V j "" ( " Vn)~ I ( \/V 

After the straightforward calculation, the error bound comes from that of (54). 


+ Yt 


dy. 


□ 


6 Vanishing Cauchy Transform 


In this section, we consider a general potential V that is defined in the sentence containing (8). 
Lemma 6.1. For a fixed z G C, let the kernel (6) satisfy the asymptotie behavior, 

\Kn{z,w)\ = ^e-'i\^-^\"(l + 0{e-^^)), 

uniformly over w G Bs{z) for some e > 0 and 6 > 0. Then there exists e' > 0 such that 

[ \Kniz,w)\'^dA{w) = 

Jc\Bsiz) 

Proof. From the reproducing property of the kernel we have 

[ \Kniz,w)\‘^dA{w) = [ \Kniz,w)\'^ dA{w) - f \Kniz,w)\'^ dA{w) 
Jc\Bsiz) Jc JBsiz) 

= Kn{z,z)-f \Kn{z,w)\‘^ dA{w). 

JBsiz) 


Taking z = w m. (59) we have 

\Kn{z,z)\ = Kn{z,z) = ^ (l + 0(e“^")). 

Integrating \Kniz,w)\ over the disk Bs{z) we have, for some ei > 0, 

[ \Kn{z, w)\^dA{w) = ^ (1 + 0{e-^^)) dA(u;) 

JBsiz) JBsiz) 

N 


TT 


(l + 0(e-^'5'))(l + 0(e-^^)). 


Proof is done by substituting the above two estimates into (60). 


(59) 


(60) 


□ 
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Lemma 6.2. For a fixed z G C, let Q be analytie in a neighborhood of z and let the kernel Kn (6) satisfy 
the asymptotie behavior, 


\-Mn{z,w) = + 

A' TT ^ ^ 

uniformly for w G Bs{z) for some e > 0 and 5 > 0. Then there exists e' > 0 such that 

dKn{z,z) = C)(e“"'^). 

Proof. We define the pre-kernel for a general Q by 


ICn{z,w) 


n—1 _ 

j=0 


(61) 


which is analytic in z and anti-analytic in w. Since ICn{z, z) = Kn{z, z) for all we have 

d 


dKn{z, z) = dlCn{z, z) = —ICn{z, w) 

OZ 

From the relation between KLn and Kn (see (18)) and from (61), we have that 

lCn{z,w) = ^(^l + 0(e~^^)^, w G Bs{z). 


Using the analyticity of ICn{z,w) in z, choosing C to be a circle around z with the radius 6/2, we have 
that 


^ICn{z,w) 


1 / K^n{C,z) 


27ri Jc (C - zfi 


dC = 0{e-^^). 


Here we use that ICniC-, z) = lCn{z, C). □ 

Proof of Theorem l.f. We recall that the probability density function for n particles is given by 

Probn({2;i, • • • ,2;n}) = \zi- 

” l</<A:<n '' j=l 

For 1 < m < n the m-point correlation function is defined as 

1 r ^ 

R/f^izi, . . . Zm) = -, - I PT0hn{{zi,--- ,Zn}) dA{zj). 

{n-my.Jcn-m ' 


We note that Rn\z) = Kn{z,z) = N pn{z) and Rn\z,w) = Kn{z, z)Kn{w,w) — \Kn{z,w)\‘^, see for 
example [22]. 
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Taking the holomorphic derivative of Pn{z) one gets 

« n - n 

dpniz) = ^ NdV{z) + '^ _ ^ Vxohn{\z,Z 2 ,--- ,Zn}) d^( 2 :j) 

J (^n 1 _ Z Zj _ 


i=2 


i=2 


= -NdV{z)pn{z) + [ ^^^Piohn {{z, W,Z 3 ,- ■■ ,Zn} Y{dA{z 


N Jc ^-1 z — w 


i=3 


Kniz,z)Kn{w,w) - \Kniz,Wp). 


In the second equality we used that our probability density function is symmetric under permutation. 
Dividing by pn we get 


dpnjz) 

Pn{z) 


= -NdV{z) + 


Kn{w,w) 1 


/C z-w 


_ f \Kniz,w)\‘^ 

Kn{z,z)Jc Z-W 


dA{w). 


(62) 


From pn{z) = jfKn{z, z), Lemma 6.2 says that the left-hand side of the above equation is exponentially 
small in N. 

Below we show that the last term of the right-hand side of (62) is also exponentially small in N. 


r \Kniz,w)\'^ 
Jc Z-W 

For the first term we find 


dA{w) = 


>Bs{z) 


\Kniz,w)\" 


Z — W 


dA(r(;) -f 


\Kniz,w)(‘‘ 


IC\Bsiz) 


z — w 


d^(ui). (63) 


'Bs{z) 


Z — W TT^ C V ^ V 


iV2 


'Bsm C 

5 p27r 


e-^e^-Nr^ 


TT- JO JO 

For the second term in (63) Lemma 6.1 says that, for some e' > 0 


(l + O(e-"^))d0dr = 


L 


C\Bs{z) 


\Kniz,w)\‘^ 


Z — W 


dA(t(;) 




'C\Bs{z) 


\Kn{z,wp dA{w) = 0{e 


We write the equation (62) replacing n by n -|- 1: 


= -NaV{.) + / Li±i(!!W iA{u.) -1- 

Pn+l{z) Jc Z-W 


\Kn+liz,w)\‘^ 


Kn+l{z,z)Jc Z-W 


d74(t(;). (64) 


By the same arguments, the left hand side and the last term in the right hand side both vanish exponen¬ 
tially fast in N. Subtracting (63) from (64), the term —NdV{z) cancels, and one obtains that 


f Kn+l{w,w) 
Ic Z-W 


dA(rc) — 


Kn{w,w) 


Z-W 


dA{w) = 


\pn{w)\^e 


2 „-NV{w) 


Ic Z-W 


-dA(tc) 


vanishes exponentially fast as N goes to infinity. 


□ 
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Discussions 


7.1 Implication of Theorem 1.1 

It is immediate that Theorem 1.1 follows from Theorem 5.9. On the other hand, a geometric intuition 
suggests that the converse is also true, i.e., the density in the normal direction from the boundary 
determines the density in any direction. In this section, we calculate the density in the tangential 
direction from the density in the normal direction and from the local shape of the curve dK. This means 
that Theorem 5.9 and Theorem 1.1 follow from each other. 

Let r be a smooth Jordan curve and 7 : M —?■ T (si—?- 7 ( 5 )) be the arclength parametrization. Defining 
n to be the normal vector 

n(7(s)) := 

let us consider the following mapping. 

H : —)> C, H(x, s) = 7 ( 5 ) + x 11 ( 7 ( 5 )). 

The Jacobian of this mapping can be evaluated to be 

1 + k(7(s))x (65) 

where k(j(s)) is the curvature of T at 7 ( 5 ). For a sufficiently small |x| the Jacobian is strictly positive 
and, therefore, H is locally invertible around 7 ( 5 ) for any s G M. It means that there exists a neighborhood 
of 7 ( 5 ) such that defines a coordinate system on that neighborhood. 

Let zo = 7 (so) £ r. We have H(0, sq) = zq- Since H is locally a diffeomorphism, there exists a neighbor¬ 
hood Uo of zq and a neighborhood Vq C of the origin such that 

Ho : Vb Uo, {x, y) ^ Ho(x, y) := H(x, sq + y) 


is a diffeomorphism. 

We define the mapping To : 1^^ —>• C by 

(T, T) ^ To(X, T) := zo + Xn(zo) + Yt{zo), t(zo) = in{zo). 

It is clear that To defines a flat coordinate around zo- 

The relation between the two coordinates, {X,Y) and {x,y), is given by the following proposition. 
Proposition 7.1. For (x,y) G Vo, there is a unique {X,Y) G such that 

~o{x,y) = Fo{X,Y). 


24 


As z goes to zq, the two coordinates (x, y) and {X, Y) are related by 


d. 


Kd, 


x = x + -Y^ + - —xy" + —x"y" - -^xy^ + ^ m - ^oH, 


dt 


K K 


y = Y - kXY + k'^X^Y - —XY^ - —^ 


d. 


- k^X^Y + 


SndsK 


X2y2 _ 


7 ndsK 
24 


y^ + - 


24 


dU 


XY^ + 0{\z-zo\^), 


where k and its derivatives are all evaluated at zq €T. 

Proof. A Taylor expansion gives 

7(y + So) = zo + i{so)y + ]p”{sQ)y^ + ^7"'(so)/ + + 0{\y\^) 

= ZQ + ty- ^{dsKii + r?t)y^ - + SAcS^Kt - K^n)/+ C>(|y|®), 

where, in the last expression, the Taylor coefficients are evaluated at zq. In the last equality we used the 
following elementary facts. 

ds^ = t, = Kt, clgt = —Kn, 

where n and t are understood to depend on 7 ( 5 ). Similarly, one can get 


n(7(y + So)) = n{zo) + 9^n(zo)y + l-d‘^n{zo)y‘^ + ^d^ii{zo)y^ + 0(1/1) 

2 0 

= n + Kty + ^{dsKt - /n)y^ + - SkOsKu - Av^t)/ + 0(|/|), 

where all the coefficients are evaluated at zo- 
Then we can expand Ho(x,y) as follows. 


^2 ^ „„,2 ^ 1 , 3 02 , \„,4 


Eo{x,y) = zo + n{zo)(x--y - ^xy - —y + --{ndsn)xy 


1 


Kd.K 


1 


+ t(2o) ( y + Atxy + ^dsKxy‘^ - —y^ - ^y"^ + -(5^At - K^)xy^ ] + 0{\z - 2o|^) 


( 66 ) 


= 2:0 + n(zo)X + t( 2 ;o)y. 

So we have obtained the series expansions of X and Y in terms of x and y. To find the series expansion 
of X and y in terms of X and Y, we simply substitute the most general Taylor expansions, 

x= ^ c,fcX^y'= + 0(|z-zo|''), y= Y 1 djkX^Y'^ + 0{\z-zof), 

j+k<4: j-\-k<4 

into ( 66 ) and determine the coefficients CjkS and dj^s by matching each order. □ 


25 






Below we illustrate how Theorem 5.9 is obtained from Theorem 1.1. We only consider ^ = 0 (i.e. purely 
tangential direction) to make the presentation simple. 

For zq = 7(so), we define x and y to satisfy the following relation. 

zo + ^ = ^o(x,y). 

Using the above x and y, we define To := 7(59 + y)- Then the identity, 

rjt 


Pn\Zo + 


y/N 


= Pn {zo + xn(To)), 


(67) 


is simply the density at the same point expressed in two different coordinate systems. This equation gives 
the density in the tangential direction expressed in terms of the density in the normal direction. This 
already implies that Theorem 1.1 gives Theorem 5.9 as we mentioned in the beginning of this section. 


Proposition 7.1 gives the asymptotic expressions of x and y in terms of p. 

dgK, V 


2 6 Ar3/2 


V = 


^/N 3 A^3/2 




Using these series expansions with 

k{zo) = k{zo) + dsK{zo) 
and the scaled normal coordinate: 


p 


up 


y/N 

dsK p^ 




2VN 6 N 


1 

iV3/2 


Theorem 1.1 gives 
pt 


2 ', 


=-2«2 


e k{zo) / ^ - 1 


1 


1 


V2tt^ y/N 




27r 'J2'K^ V 

This is exactly what Theorem 5.9 says for ^ = 0 


/ k{zq){ 1 + 2>p^) dsK{zo){p'^ + p) 
sVN 3N 




( 68 ) 


7.2 Orthogonal polynomials from density 

We have obtained the asymptotic expansion of the density using the asymptotics of the orthogonal 
polynomials. We can also do the opposite, i.e., we can obtain the subleading correction of the orthogonal 
polynomials assuming certain behavior of the density. In fact, we obtained the expression in Lemma 5.5 
by this idea. The proof in the appendix is not very motivating as it does not indicate how we find the 
identities in the first place. Below we show another proof of Lemma 5.5 that does not require any prior 
knowledge about h-i and /iq. Instead, we will use the behavior of the density given in Theorem 1.1 and 
Theorem 1.2 (or its weaker version Lemma 4.2). 
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Alternative proof of Lemma 5.5. One can integrate (58) in the tangential direction to get 


Pni Z0+ 


rjt 

Vn 


= Pn (zo) - 


1 / 1/ Im(/i_i(zo) -/io(^o)) 




+ °'iV 


1 


/ K(3r/2 + 1) 


2vr V SVN 


l(^ lm{h-i{zQ) - ho{zo)) p^ds 
\f;'izo)\ 3 




Comparing the term of order 1/N with that of (68) we obtain the second equation in Lemma 5.5. 

Let zq £ dK and zi £ IntiL be given as in the proof of Lemma 5.8. Let 0 < t < n < 1/6. Let 
Z 2 = zq + such that it is located symmetrically to zi with respect to dK. Using the same 

method as in the proof of Lemma 5.8, one can integrate (57) to get 


Pn{zi) - Pn{z 2 ) = - + ^ 
TT TTiV 




12 18k2 


dsKlm {h_i{zo) - ho{zo)) 
3\lpfzo)\K 


24k 




From Theorem 1.2 the right-hand side must be equal to l/vr up to an error given by 0(A^^/®e“'^^^). 
Taking n < 1/18, this implies that the 1/A^-term above vanishes. This provides the first equation of 
Lemma 5.5 when combined with the second equation that we already obtained. □ 


We remark that the two relations in Lemma 5.5, combined with the asymptotic behavior of hQ{z) and 
h-i{z) as z goes to oo, are enough to fully determine ho and /i_i. Lemma 5.5 gives Re(/io(2) + h-i{z)) 
and Im(/io(^) — /i_i(z)) on the ellipse. Since both functions are harmonic on the exterior of the ellipse, 
and since both are bounded at oo, Lemma 5.5 determines both functions on the exterior of the ellipse. 
This uniquely determines ho{z) + h-i{z) up to an imaginary constant and ho{z) — h-i{z) up to a real 
constant, on the exterior of the ellipse. The undetermined constants can be fixed by the behavior of 
the orthonormal polynomials’ leading coefficients pn and 7n-i (see Lemma 2.2) when n, N ^ oo while 
n = NT. 


7.3 The number of particles outside K 

In the scaling limit, most of the Coulomb particles are supported on iL. A typical sample looks like in 
Figure 3. In this section we will evaluate the number of particles outside K. This problem was suggested 
to us by Arno Kuijlaars. 

Note that the total number of particles is given by n = Ai x Aiea{K)/TT. Therefore, the expected number 
of particles outside K is given by 

^T-out ■= n — ^(particles inside K) = N [ ( -/3„(2;) ) dA( 2 ;). 

Jk / 

Inside K and away from ^/^-neighborhood (that we denote by U]\i,t in Theorem 1.2) of dK, Theorem 
1.2 says that 7r“^ — pn{z) disappears faster than any algebraic power of N in the scaling limit. This leads 
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Figure 3: Eigenvalue distribution of a sample of a random matrix Mn^t,T (71) with n = 1024, t = 0.75 
and T = 1 


to 


nout = N [ - pn(z) )d^(z) + 0{e ) 

r-L /-O 


Jo J-N'' 




1 + 


Vn 


deds + 0(e 




where L is the length of dK, 7 ; [0, L) —)• dK is the arclength parametrization of dK, and the factor 
(1 + K^/y/N)/\nV is the Jacobian coming from the change of the measure from d^ to d^ds. 

At a boundary point zq G dK, Theorem 1.1 says that 

p„(^Z0 + ^n) = ^£1(0 + + ipg'(0 + O (w-3/2+9- 


where 


= ^erfc(\/20, 


^a)(^) = -(- 0 ) ,,2 


(r - 1) e 


-252 


3^2^ 

v/^V 18 V 9k2 12k J\ 


We get 


L ^0 


^out 




= Vn 


-Af" 


rL rN‘' 


1 - p!S{() - - ^pg’(f) 


1 + k(7(s)) 


Vn 


d^ ds + o 


Vn 


>0 Jo 


p'Sio - 


(«pS’({) 




^(7('S)) 

Vn '' ^ 3V2V 

+ ^ (k(7(s))^/0^^(0 + 

V^ V [i I n ^ /^ ^^(7(g)) idsVlis))f Qf^(7(g)) ^l n 

(27r)3/2 Jq A^ \ 12 18k2(7(s)) 24k(7(s)) 


d^ ds + o 


Vn 


Vn)' 


( 69 ) 
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where we have chosen v such that 0 < lOz^ < 1/2. In the second equality, we made the substitution 
^ ^ and used the identity and symmetries: 

TT 

One notices that the subleading term of order in the penultimate formula for riout cancels out. 

And the “arclength density of escaping particles” over dK is uniform in the leading order. 

Parameterizing the ellipse by (/(e*^) with 6 G [0,27r) where the arclength element is given by ds = 
|i/)'(e*®)|d0 and using the explicit expressions for k, dgK and d^K (see ( 88 ), (89) and (90) respectively), 
the integral over the arclength in (69) can be evaluated to give 


^/n (-i + 9t"^ — 9t^ — 1 — 6 t(t^ — 1 ) cos( 20 )\ ^ ^ 

(27r)3/2 _/q Y 12 n(l + t 2 - 2tcos(26'))3 / Vl - ^\V^) 

y/Ti /4(l + t)i/2E(A:) 1 (l + t2)E(A;) + 2(l-t)2K(A;)\ /1\ 

“ (27r)3/2 (l-t)V 2 n 9(l_t)3/2(l+i)l/2 ^ 


with k 


2 y^ 

1+f 


K(k) and E(k) are the complete elliptic integrals of the first and second kind. 


K{k) = 



dO 

Y'T^^A)2~sii?^ 



1 — /c2 sin^(0)d0. 


Note that we can write the leading term also using L, the circumference of the ellipse. 


^out 


y/nL 

(27r)3/2\/r 




This expression is universally true, for any potential that shares the density profile at the boundary, 
when we define L to be the length of dK. 

The similar calculation can show that those “escaping particles” did not go far. Below we calculate the 
expected number of particles outside K in the A^^^^/s.j^gjgEi^orhood of dK: 



.( 0 ) 

roo 


(0 + 


^(7('S)) 

Vn 




g ^-1 

3\/27r3 


1 

+ iv 


(k(7(s))^pW(0 



d^ ds + o 



Comparing with the second line of (69), we notice that the “arclength density of the leaked particles” 
over dK is exactly same as that of the escaping particles since the two integrands only differs by a minus 
sign in front of the term of order 1 /^/N that turned out to be zero. 

Based on the observations in this section one may say that, when N is finite, (the normalized measure 
of) the particles are shifted from the limiting distribution 7 lie at = 00 , only around the boundary 
dK and only in the normal directions of dK. 
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re 

samples 

^out 


El 

y/n{Ei - flout) 


2 

2.41 • 

10« 

0.66073 ± 8 • 

10“^ 

0.66226 

(21.7 ± 1.2) • 10-4 

t - 0.5 

4 

1.13- 

10^ 

0.95741 ± 5 • 

10-5 

0.95822 

(16.3 ±0.9) • 10-4 


8 

7.80- 

10® 

1.36998 ± 7 • 

10-5 

1.37043 

(12.8 ± 1.9) • 10-4 


16 

2.7 • 

10^ 

1.94876±4- 

10-5 

1.94890 

(6.0 ± 1 .8) • 10-4 

r = 1 

32 

9.54- 

10® 

2.76381 ± 9 • 

10-5 

2.76382 

(0.3 ±5.0) • 10-4 


64 

1.65 • 

10® 

3.9141 ±3- 

10-^ 

3.91404 

(-4.1 ±20.2) • 10-4 


4 

2.16- 

10® 

1.3261 ±1- 

10-4 

1.33977 

(27.2 ±0.2) • 10-5 

t = 0.75 

10 

1.33- 

10® 

2.1535 ±2- 

10-4 

2.15936 

(18.4 ±0.6) • 10-® 

T = 0.4375 

32 

3.69- 

10® 

3.8944 ±2- 

10-4 

3.89639 

(11.1 ±0.9) • 10-® 


64 

2.43- 

10® 

5.5201 ±3- 

10-4 

5.52114 

(8.1 ±2.0) • 10 -® 


Table 1: Numerical simulation of random matrices flout is the number of eigenvalues outside the 

ellipse averaged over all samples, Ei the expectation value according to (70) including the first subleading 
correction. In the flout column and in the last column we show 95% confidence intervals. 


We can evaluate the number of particles outside the ellipse numerically using Matlab. We use the fact 
[25] that the particles distributed according to (2) share the same distribution as the eigenvalues of the 
random matrix ensemble, given by 

a/T 

Mn,t,T = , -^ (\/rTtAn + (71) 

V 1 — 

where and are independent copies of an n x n GUE, i.e. Hermitian matrices where the diagonal 
entries are i.i.d. AA( 0 ,l/( 2 n)) and the upper triangular entries are i.i.d. AA(0, l/(4n)) + iAA(0, l/(4n)). 
Our numerical calculation seems to support the asymptotic formula (70). We can see from Table 1 that 
the numbers in the last column tend to converge to zero as n grows (except for a few higher re’s where 
the error bounds are too large to tell). 

We also looked at the distribution of the eigenvalues outside over the arclength of the ellipse. For this we 
projected the eigenvalues outside of the ellipse onto the ellipse, in a direction normal to the boundary, 
i.e. we associated an eigenvalue outside of the ellipse with its nearest point on the ellipse. See Figure 4 
for the plot of the numerical simulation. 


A Plots of density: higher corrections 


Here we collect a few plots of pn across the boundary that supports Theorem 1.1. 

The leading behavior is shown in Figure 2. In Figure 5 we plot the subleading corrections of pn, i-e. 


pW(e) :=v/iV 

p^n\0 ■■=N 


+ -^erfc(^/20 


1 






- 2^2 
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-3-2-10123 

Figure 4: Histogram of the distribution of outside eigenvalues regarding arclength where t = 0.5, T = 1, 
n = 64, ^ of samples = 1.25 • 10®. The histogram shows the arclength density of the outside eigenvalues. 
We divide the ellipse by 140 arcs of the same length, and counted the average numbers of “escaping” 
eigenvalues that are projected to each arc. The dotted line gives the leading theoretical behavior, while 
the solid line is the theoretical curve that includes the first subleading correction. The 0 corresponds to 
the point of the maximal curvature on the ellipse. 




Figure 5: pn'^ (left) and (right) for n = 10,100,1000 and oo (dashed lines) 


According to Theorem 1.1 we know that 


lim pW(0 


lim 




e-2«^ f 2 2^5-8^3 + 3^ , 

V27r3 V 




(72) 


Both pictures are obtained for t = 0.5, T = 1 (i.e. N = n) and zo = ~ 1.73. We also mention 

that k{zo) Ri 5.19, dsK{zo) = 0 and dgK{zo) ^ —374. The gray, darker gray and black lines correspond to 
n = 10, n = 100 and n = 1000 respectively. The dashed lines are the theoretical limits shown in (72). 
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B WKB approximation 


The Hermite polynomials fulfill the differential equation 

H'l{z) — 2zH'i^{z) + 2kHk{z) = 0, A; = 0,1, 2, • • • . 


By Theorem 2.1, the orthonormal polynomial, pk, satisfies 

+ 2kpk{z) = Q, /c = 0,1,2, ••• . (73) 

We note that n = NT. 


In [10] the existence of 1/N expansion of the orthogonal polynomials, for a more general case, is shown 
using Riemann-Hilbert method (see, for instance. Theorem 7.10 in [10]; also see Section 3 in [15]). 
Using this fact, the 1/N expansion is obtained by simply positing the 1/N expansion of the orthogonal 
polynomial and plugging into the second order differential equation (73). 

According to [10, 15] we know that for I G N there exists the asymptotic expansion 

/ N \ f 111 \ 

Pn+Az) = j exp i^NTY_i{z) + Yo{z) + + j^^ 2 iz) + • • • + j^r^Yiiz) + 0{N-^-^)j 

uniformly in a compact subset of C \ [—Fqj-A^o]- (We put the prefactor so that it matches the 

behavior of 7n+r in Lemma 2.2.) Note that the functions Y^’s depend on r and T (= n/N). Putting the 
asymptotic expansion of Pn+r into the equation (73), we get 


“ = ^ 1 (*'> + E (") + E 


fc =0 


k=0 


yjM. 

ATfc 


2z(nTYU{z) + Y. 


k=0 


YM 

ATfc 


+ 2iVr + 2r + C>(iV'"^). 


Expanding in large N, one gets the following equations at each order. 


N^ : 


- 2zY/_^{z) + 2 = 0, 


N^ : ^ i2YY(z)Y/iz) + Y^iz)) - 2zY/{z) + 2r = 0, 


iV-: ^ 

2T 


(2TYU{z)Y;{z) + Y;Y{z) + Y,yiiYY',-k-Az^ -2 zY;{z) = 0, l<j<l. 

V fc=o / 


One can notice that the above gives the recursive relations to find Yj for all j > —1. The first equation 
gives 


(i - , 


(74) 
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where the sign of the square root is chosen by the asymptotic behavior: = 0{l/z). Using the 

equations at the orders and N~^ we get 


n . _ z - z^\ - (1 + 2r) 

2(z2-F2) 

4F2z2(^(l-yFnf7^)(l + 2r)+r(r-l)) + ^^(1 - 4r(r + 1)) 
~ IGTzy^l -F2/z2(z2 - F2)2 

Using Lemma 2.2, we get the following boundary behaviors. 


Y_i{z) = logz 




+ ^*(1/ z), 


Yoiz) 

Yiiz) 


-logz - (r + 2 ) logcap(Lsr) + 0{l/z), 
1 + 6r + 6r2 

--+ 0 ( 1 / 2 ;), 


Using these, one can obtain, by integrating (74) and (75), the followings. 


(75) 


y-iiz) 


Yoiz) 


Yiiz) 


^2(1 - yi - F2/^2)2 ^ J 2^(1 _ _ F,yz2) ] £ 

2F^ F(2 J 2T' 

5 (5+ 2^1-Fi/z^ ) ^^ - (’■+0 iogc«p(^o 

3Fq ^2y^1 — F^jz^ — 1 + 4r(l + — Fq 7^2 + r)^ — 222(1 + 6r(r + 1)) 

48rv/l - F2/^2(^2 _ ^2) • 


Comparing with (12) this proves Lemma 3.1. Note that Y)- depends on r if j > 0, i.e. Yj( 2 ) = Yj{z;r), 
and we have g{z) = Y-i{z) + and hr{z) = Yi{z;r). 


C Proof of three lemmas 

C.l Proof of Lemma 5.4 

From the general relations k = d arg n(s)/ds and 

n(^o) = |7^I^'(2 ;o)| (76) 

for a general smooth curve tp • ® C (see Section 5 for exaxmple), one finds that 

k{zo) = |V^^(2o)|Re ^1 - = W{zq)\ -Re(^n(zo) ^,l^°l ^, (77) 

for any boundary point zq G dK. 
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One can also find the following general relation on the boundary (we suppress the arguments): 


dsK = (in9 — in(9)|V^'| Re ( 1 — 


1p' 


12 


= Im I n 


V:"' 

■0' 


3 ir? 


2 


f\2 


For any vectors vi and V 2 the following identities hold at zq G dK. 






d-s/2 9vi 


\i^'izo)\ 

i>{zo) 


i^{zo) 

IV''(^o)| 


2Re I vi 


'ip^izo) ^ 


'4^izo) 

+ Re ( V 1 V 2 


Re vi 


2V^'(zo) 

r{zo) 


'ipizo) J 

'il^'izo)\ 

— Vi—^—- I 


'4’'{zo)J ’■'V’(^o)y 

r'izo) irizo)? 


Re V 2 


rjzo) 

i’'izo) 


V2 


-ipizo) 


'^'{zq) {i}'{zo)f 


^ ^ V (V’(^o ))2 V’(^o) } 


( 78 ) 


(79) 


In Lemma 5.4 one can see that, for any non-negative integers k and Z, 


®k,i{zo) 


V’(^o) r.kal \i^'{Zo)\ 
|V^'(zo)| " * V^(zo) 


comes from the Taylor expansion. And the uniform bound stated in the lemma comes from the real- 
analyticity of the function |V’^|/V’ around dK. Explicit calculations using (77), (78) and (79) give the 
following list of identities that hold at the boundary dK. 


© 1,0 

© 0,1 


©2,0 

© 1,1 


©0,2 


-K, 


— \'ip'\ Im 


2 fr' {r? 


+ Re nM ^ - 


ijj' {'fp'Y 


-|- \'4)'\k — i\'ip'\^ Im 


{ip'y 


- dsK + 2i\'ip'\K, 

(i” (^) + ‘) - R* - 1'*'!* + I” (I? 


All the above identities are true at the boundary. We suppress all the arguments for simplicity. 
To prove Lemma 5.4 we only need to show the following relations: 

Lemma C.l. Let zq G dK. We have the following identities 

'ip{zo)'ip"{zo) \ _ 1 dsK 

{f^'izoW ) 3\'lffzo)\K 
2 f r'jzo) _ irizo)r \ 

\i^'{zo) {fj'izoffJ 


Re n 


9 (dsK)'^ d'^K , ,,, ,, 



(80) 

(81) 
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Proof. Unlike the previous relations in this section, the above two hold only for elliptical droplets, and 
the proof can be done using the explicit parametrization of the boundary by 


7 : [0, 27r) C, 7(6*) = 0(e*®) = 
It is straightforward to calculate 


1 


e*® + te-*^ 


\/i _ t 2 (1 ^2*0)’ y 1 _ t 2 gSie ’ 

The following general relations are derived from 'il>(4>(u)) = u. 

Using (76) and the above two sets (83) (84) of equations one gets 

pie _ fp-i9 

n(7(0)) = = 

+ ^2 — 2t cos(20) 


i'll f„i9 


in = - 


T 6t 
1 — e^*® 


(/)'(u)^ (j)'{u)^ 


Using (83) (84) (85), one obtains the following expressions. 


n{-f{e)f 


m^)Wiin) 

i’'iii0)) 

riije)) irii{e))? \ 

^'iim if^'iimf) 


2t {t - e-n 

| e *® — te “*^|2 ’ 

2t(l -t2)(3 + te-2i^)(e-® 
T |e® — te“*®|® 


te*^)2 


(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 


that appear in the left-hand sides of (80) and of (81). Taking the imaginary and the real part respectively, 
one gets the following evaluation at 7(0) G dK. 


Re 




(V’^^y / 


7(0) 

7(0) 


2t sin(20) 

1 -|- t2 — 2t cos(20) ’ 

2t(l — f2)(i3 _ _|_ ^2 cQg(|26*) -|- 3cos(20) -|- fcos(40)) 

r|l -bt2 - 2tcos{29)\^ 


The proof will complete if we can show the following relations: 

/f _^2\3/2 

nm = ^ ’ 


\/T(l + t2 

— 2t cos(20))^'^^ 

6(1 — sin(20) 

T(l + t2 

- 2tcos(26l))^’ 

12(1 -t2 

)®/^t(2t(cos(40) — 


r3/2 (1 +t2 _ 2t, 


1 


y/l + t'^ — 2t cos(20) 


( 88 ) 

(89) 

(90) 

(91) 
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that appear in the right-hand sides of (80) and (81). We remind that the function 7 is defined at (82). 

The first relation ( 88 ) is obtained using ( 86 ) and (77). The second (89) and the third (90) are obtained 
by the direct calculation using 


dsK{-f{e)) 


1 dK(7(6>)) 

|(/)'(e*®)| d9 




1 d{dsK{j{e))) 

|0'(e*^)| d0 


The equation (91) is simple to obtain using (83) (84). 


□ 


C.2 Proof of Lemma 5.5 


Using (25) one obtains that 

hohi^)) + h-i{'y{e)) = - 


3F^ + 2z^ 


2ATy/l - F^/z^ (z2 _ ir2) 
Taking the real part, this gives 


(e*® -I- te (e^*® -|- t^e -|- 8t) 


z='y{e) 


12T -te 


.id _ -l-a — id^^ 


Re (hoi'jid)) + /i_i(7(6'))) = - 


(1 — (1 — 26t^ + t'^ + 6 f cos(20) -|- 6t^ cos(26l) -|- 12f^ cos(40)) 

12r(l_ 2tcos(26'))^ 


(92) 


which is the left-hand side of the first equation to prove. Using ( 88 ) (89) one can confirm the first equation 
of Lemma 5.5. 


Using (25) one obtains that 

h-i{-f{e)) - ho{-/{0)) = - 




4T(z2 - 


t (e-2*e + fe^^o - 2t) 


2=7(0) 


T le*® — te 


Using (88)(89)(91) one can confirm the second equation of Lemma 5.5. 


C.3 Proof of Lemma 5.6 


Using (85)(88)(89)(91) we get 


VT{nt - n) \i^'{zo)\ 
Vl - F i^izo) 
VT{m + n) \i^'izo)\ 
Vl - F V'(^o) 


l + t‘^ - 2te-2*^ _ ^ idsK 

1 + F — 2tcos{29) 3\^jJ'{zo)\K' 

1 — K 

1 + F — 2t cos{29) \ip'{zo)\ 
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